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1.  Introduction 

Understanding  of  localized  states  on  the  materials’  surfaces  is  crucial  to  understand  a 
multitude  of  physical  and  chemical  phenomena  occurring  on  the  surfaces.  Owing  to  the 
pioneering  works  by  Tamm  and  Shockley  and  following  numerous  efforts,  the  origin  and 
nature  of  the  surface  states  has  been  revealed  [1,2].  Optical  analogues  of  the  surface  states 
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also  have  been  studied  in  photonic  lattice  systems  such  as  photonic  crystals  [3-6]  and  optical 
waveguide  arrays  [7,8].  The  photonic  lattices  provide  ideal  and  controlled  environments  for 
investigating  surface  modes,  which  have  enabled  new  approaches  to  photon  manipulation. 
Localized  modes  also  can  be  supported  on  the  surface  of  metamaterials  made  of  metal- 
dielectric  multilayers,  leading  to  deep  subwavelength  mode  profiles  [9].  Such  surface  modes 
supported  either  in  solids  or  in  photonic  lattices  are  formed  by  termination  of  the  periodicity 
and  can  be  categorized  into  two  types,  Tamm  and  Shockley  states:  while  Tamm  states  appear 
due  to  strong  perturbation  above  a  certain  threshold  at  the  surface,  Shockley  states  requires  a 
different  condition,  break  of  the  stronger  bonding  (coupling)  at  the  surface  in  the  chains  of 
two  alternating  weak  and  strong  bondings.  In  a  band  diagram,  Shockley  states  show  up  in  an 
inverted  bandgap  which  is  formed  with  crossing  of  adjacent  bands. 

The  metal-dielectric  metamaterials  offer  intriguing  platforms  for  investigating  surface 
modes.  Weak  and  strong  (binary)  couplings  for  formation  of  Shockley-like  modes  can  be 
realized  not  only  with  different  strengths  but  also  with  combinations  of  signs  for  the  two 
coupling  coefficients,  which  have  not  been  realized  in  the  previous  configurations  [9].  In  this 
paper,  we  present  the  origins  of  the  surface  mode  formation  arising  from  termination  of  the 
periodicity  and  identify  their  existence  conditions  in  metal-dielectric  metamaterials.  We 
describe  how  internal  interactions  and  surface  perturbations  affect  the  formation  of  surface 
modes  in  semi-infinite  metamaterials  through  a  simple  model,  which  also  allows  clear 
physical  insight  into  the  localized  modes  in  general  photonic  systems. 

2.  Internal  interactions  &  band  crossing  in  binary  metamaterials 

We  begin  by  describing  unique  bulk  properties  of  infinite  binary  metamaterials  in  Fig.  1(b) 
and  1(c)  via  comparison  with  conventional  metal-dielectric  metamaterials  in  Fig.  1(a).  The 
metal-dielectric  multilayers  in  Fig.  1(a)  have  been  studied  as  metamaterials  for  sub  wavelength 
imaging  and  focusing  since  we  can  tailor  their  dispersions  at  will  by  adjusting  layer 
parameters  [10-14].  For  the  conventional  structure  in  Fig.  1(a),  the  effective  medium 
dispersion  curve  k2  / £{l+k2  / e±  =  k2o  ,  where  s±~'  =  '  and  for  Metric  fields 

acting  along  the  direction  normal  and  parallel  to  the  layer  plane  and  f. is  a  filling  ratio  of 
the  zth  layer,  agrees  well  with  the  accurate  dispersion  calculated  by  the  transfer  matrix  method 
[15]  for  small  k  (  =  kx A ,  A  is  the  lattice  period).  However,  this  approach  can  break  down  for 

the  structures  that  have  multiple  interactions  between  constituents  in  metamaterials.  Since  the 
effective  medium  parameters  are  obtained  by  macroscopically  averaging  out  electromagnetic 
properties  of  each  constituent,  internal  microscopic  features  are  not  fully  taken  into  account. 
Figure  2(d)-2(h)  display  the  microscopic  details  for  the  coupling  between  surface  plasmon 
polariton  (SPP)  modes  in  the  layers  [16].  The  single  gap  SPP  mode  (TM  polarization)  in 
Fig.  2(d)  couples  to  the  adjacent  identical  mode  through  a  thin  metal  layer  as  shown  in 
Fig.  2(e),  which  causes  mode  splitting  into  symmetric  and  antisymmetric  mode.  For  this  case 
the  propagation  constant  of  the  antisymmetric  mode  ( pa )  is  bigger  than  that  of  the  symmetric 

mode  ( ps )  unlike  conventional  index-guided  waveguides.  Thus  the  corresponding  coupling 
coefficient  defined  by  c  =  (/?-/?)/  2  is  negative  (anomalous  coupling)  [17].  Now  let  us 

consider  the  structure  in  Fig.  1(b).  The  multilayer  is  similar  to  the  one  in  Fig.  1(a)  except  for 
the  alternating  metal  layer  thicknesses,  40nm  and  20nm,  which  leads  to  two  different  strengths 
of  the  couplings  between  adjacent  SPP  modes.  As  a  result,  two  split  dispersion  curves  are 
obtained.  We  can  easily  notice  that  the  dispersion  by  effective  medium  approximation  fails, 
even  for  small  k  .  Similar  effect  is  observed  with  the  structure  in 
Fig.  1(c),  whose  unit  cell  consists  of  one  metal  layer  (20nm)  and  three  dielectric  layers, 
consecutively  n  =  1.34  (50nm),  n  =  3.48(200nm),  and  n  =  1.34(50nm)  supporting  single 
hybrid  SPP  mode  as  seen  Fig.  2(f)  [18].  Although  the  structure  looks  more  complicated  than 
the  one  in  Fig.  1(b),  there  are  also  two  types  of  couplings  between  SPP  modes  through  the 
high  index  dielectric  layer  [normal  coupling,  c  >  0 ,  Fig.  1(g)]  and  through  the  metal  layer 
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[anomalous  coupling,  c  <  0 ,  Fig.  1(h)].  Thus,  the  resulting  two  dispersion  curves  with  huge 
splitting  cannot  be  approximated  by  the  single  effective  medium  curve  in  Fig.  1(c).  The 
examples  in  Fig.  1(b)  and  1(c)  clearly  demonstrate  the  effect  of  the  multiple  interactions  in 
metamaterials,  which  also  plays  an  essential  role  in  supporting  surface  modes  in  metal- 
dielectric  metamaterials. 
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metal 


H  dielectric 


Fig.  1.  (a)-(c)  Three  types  of  metal-dielectric  metamaterials  and  the  corresponding  dispersion 
curves  for  the  unit  cell  parameters  (a)  metal  (Au,  20nm,  dark  grey)  and  dielectric  (n  =  1.5, 
lOOnm,  light  orange),  (b)  metal  1  (Au,  20nm),  metal  2  (Au,  40nm),  and  dielectric  (n  =  1.5, 
lOOnm),  (c)  metal  (Au,  20nm),  dielectric  1  (n  =  1.34,  50nm,  yellow),  dielectric  2  (n  =  3.48, 

200nm,  light  grey),  and  dielectric  3  (n  =  1.34,  50nm).  L:  low  index,  H:  high  index  dielectric. 
k  is  the  transverse  wave  vector  multiplied  by  the  lattice  period  a  .  The  excitation  wavelength  is 
1550nm.  (d)-(e)  Normal  and  anomalous  coupling  between  two  coupled  SPP  modes  in  metal- 
dielectric  multilayers.  Blue  curves  indicate  transverse  electric  fields  ( Ex ).  (d)  single  SPP  mode 
in  a  metal-dielectric-metal  layer,  (e)  symmetric  and  antisymmetric  mode  by  anomalous 
coupling  of  two  gap  SPP  modes  through  a  thin  metal  layer,  (f)  single  hybrid  SPP  mode  in  a 
metal-three  dielectric  layer  (low  index/high  index/low  index),  (g-h)  symmetric  and 
antisymmetric  hybrid  SPP  modes  by  normal  coupling  through  a  high  index  dielectric  layer  and 
by  anomalous  coupling  through  a  thin  metal  layer,  respectively.  The  diagrams  in  the  right 
display  schematic  mode  splitting  in  the  propagation  constant. 

For  the  binary  metamaterials  in  Fig.  1(b)  and  1(c),  the  band  structure  evolves  with  the 
relative  strength  of  the  couplings  between  SPP  modes.  We  define  asymmetry  constant  as  the 
ratio  of  coupling  coefficients,  tj  =  C  /  C  .  The  asymmetry  constant  77  can  take  both  negative 
and  positive  values  depending  on  the  sign  of  C  and  C  .  For  example,  77  is  positive  for  the 
structure  in  Fig.  1(b)  and  negative  for  the  one  in  Fig.  1(c)  [9].  The  evolutions  of  bands  are 
shown  in  Fig.  2(a)-2(c)  for  the  structure  in  Fig.  1(c)  with  increasing  magnitude  of  77 .  To 

control  the  magnitude  (through  c  ),  we  adjust  the  thickness  of  the  metal  layer  to  15nm,  8.4nm, 
and  5nm,  respectively.  The  thickness  of  the  other  layers  is  fixed  (see  the  figure  caption).  The 
thinner  the  metal  layer,  the  stronger  the  coupling  becomes,  thus  \c  \  increases.  The  evolution 
of  two-band  structure  as  a  function  of  77  is  simplified  in  Fig.  2(d)  by  coupled  mode  theory 
[2,9].  With77  =  0,  we  obtain  only  two  discrete  values  corresponding  to  the  symmetric  and 
antisymmetric  modes  from  the  normal  coupling  ( c+ )  with  no  band  formation  [two  red  circles 
in  (d)].  As  the  magnitude  of  77  increases,  the  bandwidth  linearly  increases  until  ]t7|  =  1 .  The 
bandgap  with  \n\<l  is  a  direct  bandgap  [the  yellow  region  D  in  (d)]. 
Figure  2(a)  corresponds  to  this  case.  The  coupling  through  the  metal  layer  (dAu  =  15nm)  is  not 
strong  enough  yet,  so  77  is  in  the  range  of-1  <  77  <  0  and  the  bandgap  is  direct.  The  eigenmodes 
at  the  bandedges  B  and  T  are  shown  in  Fig.  2(e)  and  2(f),  respectively.  With  reduced  metal 
thickness  (8.4nm  from  15nm),  C  is  increased  and  reaches  the  same  magnitude  as  C  ( 77  =  -1 ). 
The  bandgap  is  closed  [Fig.  2(b)]  at  k  =  0  and  the  eigenmodes  in  (e)  and  (f)  become 
degenerate.  The  touching  point  becomes  a  singularity  [green  arrow  in  (b)]  called  diabolical 
point  or  Dirac  point.  It  brings  important  consequences  for  the  light  propagation  dynamics  in 
the  nano- scale  as  the  intersection  of  energy  bands  (conical  intersection)  has  great  importance 
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for  molecular  dynamics  and  electron  transport  [14].  As  77  increases  more  (|t7|>i),  the  two 

bands  cross  and  an  inverted  bandgap  appears  [grey  region  in  (d)].  As  a  consequence  of  the 
band  crossing,  the  eigenmodes  at  the  bandedges  are  interchanged. 
Figure  2(c)  shows  a  crossed  band  with  r/<- 1  when  the  metal  layer  thickness  is  more  reduced 

to  5nm  and  |c_|  further  increases.  The  crossed  band  can  be  easily  identified  by  the 

eigenmodes  at  the  bandedges.  We  can  observe  the  exchanged  eigenmode  symmetry  at  B’  and 
T’  in  Fig.  2(g)  and  2(h).  The  band  crossing  feature  in  the  binary  metamaterial  is  analogous  to 
the  electronic  band  crossing  for  the  chain  of  atoms  with  s-  and  p-  orbitals  [2].  The  band 
crossing  phenomena  is  one  of  the  essential  features  of  binary  metamaterials  and  causes  a 
direct  consequence  on  the  formation  of  surface  modes. 


Fig.  2.  Bandgap  closing  &  band  crossing,  (a),  (b),  (c)  Band  diagrams  for  the  metal  (Au)  layer 
thickness  dAu  =  15nm  (a),  for  dAu  =  8.4nm  (b),  and  for  dAu  =  5nm  (c).  The  thickness  of  the  other 
layers  is  kept  as  dMgF2  =  50nm  and  dSi  =  200nm.  B,  B’:  band  edges  of  the  upper  bands,  T,  T’: 
band  edges  of  the  lower  bands  in  (a)  and  (c).  The  arrow  in  (b)  indicates  the  singular  point,  (d) 

Band  evolution  calculated  by  coupled  mode  theory  as  a  function  of  77  .  The  bands  are  marked 
in  blue,  the  direct  bandgap  (D)  in  yellow,  and  the  inverted  bandgap  (I)  is  in  grey.  The  arrow  in 
(d)  marks  the  bandgap  closing  point  C  at  77  =  -1 .  (e-h)  Eigenmodes  at  the  bandedges  B  (e),  T 
(f),  B’  (g),  and  T’  (h),  respectively. 


3.  Formation  of  surface  modes  and  their  existence  condition 

When  the  periodicity  is  terminated,  localized  modes  whose  amplitudes  decay  away  from  the 
surface  can  be  supported.  To  analyze  the  existence  conditions  of  the  surface  modes,  we 
describe  their  formations  based  on  coupled  mode  theory.  For  a  semi-infinite  metamaterial 
with  termination  of  C  coupling  on  one  side  as  seen  at  the  top  in  Fig.  3(c),  the  governing 
coupled  differential  equations  for  the  bulk  are  given  by  Eq.  (1): 

I  i  +  j 3  j  $  4-  C  b  l  +  C  b  —  0,  f  i  +  J30  —  0,  (1) 

\  dz  )  n  "  ""  +  "  \  dz  )  n 

where  fio  is  the  propagation  constant  of  the  SPP  mode  in  an  isolated  unit  cell,  and  an  and 
bn  are  the  field  amplitudes  of  the  nth  odd-numbered  and  even-numbered  SPP  modes  from  the 
surface  of  binary  metamaterials,  respectively.  Next  we  apply  the  boundary  condition  on  the 
surface  (id  /  dz  +  p  )ax  +  Cbx  =  0 ,  where  p  is  the  propagation  constant  of  the  surface  mode,  and 
ax  and  b  are  the  amplitudes  for  the  first  and  second  SPP  modes  from  the  surface, 
respectively.  For  trial  solutions,  we  use  an  =  peipz  and  bn  =  qe'Pz  where  pn  =  prn~x  and 
q  =  qorn  l .  For  the  localized  modes  decaying  away  from  the  surface,  the  condition  |r|  <  l  must 
be  met.  By  inserting  the  ansatz  into  Eq.  (1),  we  obtain  as  below: 
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q„  =  (K  +  Z)  p0  -Kp^+^f  +  lf  =0,  -  Kqa  +{rjr  + 1)  p0  =  0, 


(2) 


where  we  introduce  new  parameters,  dimensionless  propagation  constant  and  surface 
perturbation,  K  =  (j3 -  j30)  /  C+  and  z  =  (j3o-  j3  )  /  C  ,  respectively.  Depending  on  the  existence  of 
surface  perturbation,  we  have  two  scenarios  for  the  surface  modes.  Here  we  will  describe  only 
negative  77  case  and  extension  to  positive  77  case  is  straightforward. 


1 

08 


(a) 


(b) 


(c) 


Fig.  3.  (a)  Existence  conditions  for  surface  modes  in  r,-z  plane  corresponding  to  the  plus 
solution  in  Eq.  (3).  The  contours  with  constant  r  indicate  degree  of  localization,  (b)  Existence 
conditions  corresponding  to  the  minus  solution  in  Eq.  (3).  (c)  Top:  A  schematic  for  semi¬ 
infinite  binary  metamaterials  with  termination  of  C-  coupling  at  the  left  end.  The  an  and  bn 
indicate  amplitudes  of  SPP  modes  at  nth  even-numbered  and  odd-numbered  SPP  modes  from 
the  surface,  not  the  physical  layers.  The  arrow  and  Z  indicate  the  location  of  surface  and  the 
surface  perturbation,  respectively.  Bottom:  Band  diagram  for  positive  (red  dotted)  and  negative 
(solid)  7/  ,  and  shift  of  the  surface  mode  locations  in  the  band  diagram  with  increasing  surface 
perturbation  Z.  (d)  The  corresponding  eigenvectors  for  the  surface  modes  in  the  regions  a-f 
marked  in  (a)  and  (b). 


(1) .  No  surface  perturbation  (  z  =  0  ):  We  obtain  K  =  0  ,  qo  =  0  ,  and  r  =  -I/77  from  the 

Eq.  (2).  Thus  we  have  localized  surface  modes  at  the  center  of  the  bandgap  ( p  =  j30 ) 
under  the  condition  of  I77I  >  l .  This  localized  mode  corresponds  to  the  Shockley  states 
in  condensed  matter.  This  mode  appears  only  in  the  inverted  bandgap,  meaning  that 
the  stronger  coupling  ( c  )  is  terminated  at  the  surface.  Only  the  odd-  numbered  SPP 
modes  ( a  )  are  non-zero. 

(2) .  With  surface  perturbation  (z  *  0):  For  this  case,  the  following  solution  is  obtained 

for  r  from  Eq.  (2). 


r  = 


(3) 


For  localized  mode  solutions,  the  allowed  ranges  of  two  parameters,  77  and  z  are  identified  in 
Fig.  3(a)  and  3(b)  from  the  two  solutions  in  Eq.  (3).  The  contours  for  constant  r  indicate  the 
degree  of  localization;  the  smaller  values,  the  stronger  localization.  Depending  on  the 
combination  of  77  andz ,  the  whole  plane  is  divided  into  several  regions.  The  existence  and 
character  of  the  surface  modes  is  determined  by  which  region  the  modes  belong  to.  We 
describe  the  detailed  properties  of  the  modes  in  each  region  below. 
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(i)  -1  <  77  <  0  &  z  <  \rj\  :  In  this  region  no  surface  mode  exists.  The  bandgap  is  direct  and 

the  surface  perturbation  is  too  weak  to  induce  surface  modes. 

(ii)  -1  <  77  <  0  &  z  >  \rj\  (the  region  e  and  f):  Tamm-like  surface  modes  are  supported  due 
to  strong  surface  perturbations.  Two  different  types  of  solutions  are  obtained  in  this 
region.  One  of  these  is  located  in  the  center  bandgap.  These  Tamm-like  modes  move 
toward  the  center  from  the  bandedges  with  increasing  surface  perturbations.  The 
other  type  of  solution  is  located  out  of  the  bands  (not  between)  and  move  away  from 
the  band  edges  with  increasing  perturbation.  Since  those  two  types  are  located  in 
different  positions  in  the  band  diagram,  their  mode  symmetries  are  also  different.  The 
mode  profiles  for  those  located  in  the  center  bandgap  are  shown  in  Fig.  3(d)  for 
opposite  signs  of  Z.  The  mode  energy  is  more  localized  in  the  even-numbered  SPP 
modes  (b  ). 

(iii)  tj  <  -1  &  z<\rj\  (the  region  a  and  b):  Although  surface  perturbation  is  weak, 
Shockley-like  modes  are  supported  due  to  the  break  of  stronger  coupling.  The 
increase  of  surface  perturbation  |z|  makes  the  surface  modes  more  delocalized  until 
it  reaches  the  line  of  z  =  +77  where  the  modes  become  extended,  which  is  in  sharp 
contrast  to  the  region  -l  <  77  <0  where  the  modes  need  strong  surface  perturbation 
( |z|  >  1 77 1 )  to  be  localized.  The  shift  of  mode  locations  in  Fig.  3(c)  also  confirms  this. 
The  modes  move  toward  the  bandedges  from  the  bandgap  center  with  increasing 
perturbation.  The  corresponding  eigenmode  profiles  are  also  shown  in  Fig.  3(d).  The 
mode  energy  is  localized  only  in  a  SPP  modes  without  surface  perturbation  and  the 
field  amplitude  decays  into  the  bulk  without  changing  the  sign.  With  increasing 
perturbation,  the  mode  energy  in  b  SPP  modes  starts  to  increase. 

(iv)  77 < -1  &  z>\rj\  (the  region  c  and  d):  Shockely-like  modes  do  not  exist  in  this  strong 
perturbation  region  but  Tamm-like  modes  can  be  supported.  The  modes  are  located 
out  of  the  bands  and  move  away  from  the  bandedges  with  increasing  perturbation. 
The  corresponding  eigenmodes  shown  in  Fig.  3(d)  change  the  sign  of  the  field 
amplitude  at  every  other  SPP  modes. 

The  mode  symmetries  for  positive  77  are  opposite  compared  to  the  case  of  negative  77  due 
to  the  flipped  band  curvatures.  For  finite  metal-dielectric  metamaterials  with  strong  internal 
interactions,  surfaces  at  each  side  can  interact,  thus  can  complicate  the  analysis  [9,19]. 
Nevertheless  the  simple  approach  described  in  this  paper  sheds  clear  insight  on  the 
fundamental  physics  of  the  surface  modes  in  metal-dielectric  metamaterials. 

5.  Conclusion 

We  have  analyzed  how  internal  interactions  and  surface  perturbations  affect  the  formation  of 
surface  modes  in  binary  metal-dielectric  metamaterials.  A  simple  model  has  been  introduced 
to  analyze  the  existence  conditions  and  mode  symmetries  of  the  surface  modes  in  semi¬ 
infinite  metamaterials.  The  surface  modes  in  metamaterials  could  be  utilized  for  manipulating 
surface  plasmon  polaritons  in  the  nano-scale. 
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